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The study of stable accretion disc oscillations relevant to black hole candidate (BHC) systems 
dates back over twenty years. Prior work has focused on both unstable and (potentially) stable 
disc oscillations. The former has often been suspected of being the underlying cause for the ob- 
served broad-band variability in BHC, whereas the latter has had little observational motivation 
until quite recently. In this article, we review both the observations and theory of (predomi- 
nantly) stable oscillations in BHC systems. We discuss how variability, both broad-band and 
quasi-periodic, is characterized in BHC. We review previous claims of low frequency features in 
BHC, and we discuss the recent observational evidence for stable, high frequency oscillations 
in so-called 'galactic microquasars'. As a potential explanation for the latter observations, we 
concentrate on a class of theories- with a rich history of study- that we call 'diskoseismology'. 
We also discuss other recent alternative theories, namely Lense-Thirring precession of tilted 
rings near the disc inner edge. We discuss the advantages and disadvantages of each of these 
j theories, and discuss possible future directions for study. 



1. Overview of Black Hole States 

Galactic X-ray sources are typically identified as black hole candidates (BHC) if they 
have measured mass functions indicating a compact object with M ^ 3 M@, or if their 
high energy spectra (~ 1 keV — 10 MeV) and temporal variability (~ 1CU 3 — 10 2 Hz) 
are similar to other BHC. A review of the general observations, a number of theoret- 
ical models, plus individual descriptions of approximately twenty galactic BHC can 



be found in Tanaka & Lewin (1995). A review of the timing analyses can be found 



in van der Klis (1994) and van dcr Klis (1995). (These latter reviews attempt to link 



the energy spectra and timing observations, and they draw analogies to similar observa- 
tions of neutron stars in low mass X-ray binaries.) A review of a subset of BHC with 
reasonably well-determined mass functions can be found in |Nowak (1995)1 



The energy spectra of BHC have been historically labelled based upon observations 
of the soft X-ray band (~ 2 — 10 keV). Intense, quasi-thermal flux is referred to as 
the "high" state. Non-thermal flux in this band, typically a power law with a photon 
index]]] of ~ 1.7, indicates that the BHC is in the "off" state (for extremely low intensity 
flux) or "low" state (for moderate intensity flux). The high state tends to have little 
variability, with a root mean square (rms) variability (cf. §|2|) of a few percent, whereas 
the low state tends to have an rms variability of several tens of percent. If a black 
hole candidate has a quasi-thermal soft X-ray component and significant high energy 
emission — occasionally modeled as a power law with a photon index ~ 2.5 — it is said 
to be in the "very high" state [cf. [Miyamoto, et al. (1991) | . The rms variability of the 



very high state tends to be greater than that of the high state, but less than that of the 



low state [cf. Miyamoto, et al. (1992) , Miyamoto, et al. (1993) |. Note that the labels 



"off" , "low" , "high" , and "very high" are purely qualitative in nature, but historically 
have been in popular use. No universally agreed upon quantitative definition for these 

f Photon index here and throughout shall refer to the photon count rate, such that photon 
index T implies # photons/keV/s/cm 2 cx E~ r , where E is the photon energy. 
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Figure 1. Left: The ratio of luminosity in a power law tail to total luminosity (I/nt/LtdO 
versus the ratio of total luminosity to Eddington luminosity (L/L-Edd) for black holes with 
well-determined mass functions. (Note: error bars are not plotted.) Right: Hypothe tical evolu- 
tionary tracks for changes in BHC luminosity. [Figure adapted from Nowak (1995).] 



states exist. More recently, the terms "hard state" and "soft state" have begun to be 
used in place of "low state" and "high/very high state", respectively. 

Figure [l] depicts this state behavior for several BHC with reasonably well-determined 



mass functions [cf. Nowak (1995) |. This figure shows the fraction of total luminosity 



in a power law- like tail as a function of the BHC's fractional Eddington luminosity. 
BHC tend to be dominated by a power law tail for luminosities below ~ 10% LEdd, 
be completely dominated by a quasi-thermal, disc-like component near ~ 10% LEdd, 
and then begin to show a power law-like tail for L ^ 10% LEdd- Above 10% LEdd, 
the fraction of total luminosity in this power law tail (which typically has photon index 
r ~ 2.5) tends to increase with fractional Eddington luminosity. There are, of course, 
exceptions to the trends shown in Fig. [l]; however, some transient BHC, such as Nova 
Muscae, have followed this trend very closely [cf. Miyamoto, et al. (1994)| . Historically, 
the soft flux component has been associated with a "classical" accretion disc, as discussed 



by Bhakura fc Sunyaev (1973) 



2. Measuring Variability 

The states of BHC have been further distinguished by their variability behavior, 
as discussed by van der Klis (1994), van der Klis (1995), [Miyamoto, et al. (1992) , and 
Miyamoto, et al. (1993). As shown in Figure |2j, the greater the fraction of luminosity 
in a hard tail, the greater the amplitude of the observed variability. Miyamoto, et al. 
(1992) and Miyamoto, et al. (1993) further point out that states appear to have "canon- 
ical" variability behaviour not only in terms of amplitude, but also in terms of frequency 
distribution. This frequency distribution is typically described via Fourier transfrom 
techniques [cf. van der Klis (1994)], which we briefly describe below. 

The key assumption in applying Fourier transform techniques is the assumption that 
the lightcurve represents a statistically stationary process [cf. Davenport & Root (1987)|. 
Effectively, this means that one assumes that individual data segements from the lightcurve 
are a good representation of an ensemble of statistically identical lightcurves. It is from 
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Figure 2. RMS variability versus ratio of luminosity in a power law tail to total luminosity 
(cf. Fig. hi). (Consistent energy bandpasses were not used for all depicted data points, although 
they predominantly r epresent a ~ 2 — 35 keV bandpass. Again, error bars are not plotted.) 
[Figure adapted from Nowak (1995).] 



these individual lightcurve segments that one forms the power spectral density (PSD) 
[cf. |van der Klis (1994| . 

The PSD is calculated by dividing the lightcurve into segments of equal length and 
then taking a Fast Fourier Transform (FFT, i.e. discrete transform, cf. Press, et al. 
1992) of each data segment. The squared amplitude of each individual FFT is then 
averaged together. Often one also averages over Fourier frequency bins. (Typically one 
chooses a scheme such that the number of frequency bins averaged over is oc /, where f 
is the Fourier frequency.) This yields the resulting PSD for each lightcurve. 

There are a number of different possible normalizations for the PSD; however, in 



this work we shall only refer to the normalization of Miyamoto, et al. (1992), For this 



normalization, the integral of the PSD over positive Fourier fre quency yields the square o f 
the total root mean square (rms) amplitude of the variability | Miyamoto, et al. (1992) |. 
That is, the PSD integrated o ver positive frequency yi e lds ((x 2 ) — (x) 2 ) / (x) 2 1 w here 



x 



is the photon count rate [cf. Miyamoto, et al. (1992) , Miyamoto, ct al. (1993)| . The 



rms can be calculated over specific Fourier frequency and energy spectral bandpasses. 
For a given narrow Fourier frequency interval, the rms is (to within a factor of -\/2) the 
fractional amount by which the lightcurve is sinusoidally modulated in that given Fourier 
frequency interval. 

It is this measure of rms variability, typically calculated over frequency intervals of 
~ 10 -3 — 10 2 Hz in BHC, that leads to the characterization of power law tail-dominated 
sources (i.e. low and very high state sources) as being the most variable BHC. As a 
function of Fourier frequency, the PSD of all states can usually be described by (possibly 
multiple) broken power laws. However, occasionally superimposed on top of this broad 
band variability are narrow — compared to the the broad band PSD — features which 
are referred to as quasi-periodic oscillations (QPO). The properties of these features are 
the focus of the next section. 
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3. Quasi-periodic Variability observed in BHC 

A quasi-periodic oscillation is usually taken to be any feature in the PSD that is well-fit 
by a Lorentzian of the form: 

PSD(f) = — r R2 Q fQPO T . (3.1) 



H/qpo-/) 2 + ff 



QPO 



As written in eq. 3.1, this functional fit is relevant to the normalization of Miyamoto, et 
al. (1992). In the above, R is the fractional rms, Jqpo is the QPO frequency, and Q is 
the mode 'quality factor'. If Afqpo is the full-width half maximum of the QPO feature, 
then Q s» fqpol 'A/qpo- 

For Q 3, it is difficult to attribute the QPO to a discrete feature, as opposed to 
being merely an extension of the broad-band variability. For Q ^ 10, the QPO is usually 
quite apparent in the PSD, and may be attributable to a discrete process in the BHC 
system. The question then arises of whether or not such a QPO feature represents a 
quasi-stable mode in an accretion disc. If the QPO feature is attributable to a mode or 
set of modes, then t here are many possi ble mechanisms for generating a finite feature 
width in the PSD [cf. |van der Klis (1994)| . Among these possibilities are: multiple modes 
distributed over frequencies Jqpo ± AJqpo, a driven-damped mode with damping time 
scale td ~ Q/qpo' a m °de that executes a random walk in phase of ~ 27r on time scales 
of order td, etc. Rarely does one have the statistics to determine which of these processes 
is responsible for the Q of the QPO. 

Prior to the launch of the Rossi X-ray Timing Explorer (RXTE), observations of QPO 
in BHC were made with Ginga, BATSE, and Granat/SIGMA. All of these observed 
QPO were at frequencies ^ 10 Hz, although the above instruments were for the most 
part limited to 60 Hz. In addition, the rms variability of these features were all 10%, 
with Q ^ 0(10). Below, we review some of the historical observations of black hole QPO. 

3.1. Historical Observations 



Reviews of some of the observed QPO features in BHC can be found in Nowak (1995) 
and Tanaka fc Lewin (1995)[ Here we consider separately low state and high/very high 



state observations of BHC 

In the low state, Cygnus X-l has shown QPO with rms as high as 15% [cf. Kouveliotou, 
et al. (1992a), Ubertini, et al. (1994), Vikhlinin, et al. (1994)]. The features were seen 
to have frequencies ~ 0.04 and 0.07 Hz. However, they appeared more as broad peaks in 
the otherwise broken power law PSD, and did not appear as discrete, narrow features. 
Similar features, with frequencies 0.04 and 0.2 Hz were seen in the hard state of the X-ray 
transient GRO J0422+32 (a.k.a. Nova Persei) [Kouveliotou, et al. (1992b), Sunyaev, et 
al. (1992), Sunyaev, et al. (1993)]. 

The source GX 339-4, which has been observed to transit through soft and hard X-ray 



states, has shown a relatively narrow [Q ~ O(10)] QPO at 0.8 Hz [Grebenev, et al. (1991)] 



during its hard state. On the other hand, in its very high state GX 339-4 has also shown 



a narrow QPO at 6 Hz [\4iyamoto, et al. (1991)]. Similar to this feature, the X-ray 



transient Nova Muscae has shown QPO in the range 3 — 8 Hz during its very high state 
[ Kitamoto, et al. (1992)| . 



In the soft state of LMC X-l (which is the only state that has been observed in this 
source), there has been a claim of a 0.08 Hz feature with rms ~ 4% [Ebisawa, Mitsuda, 
& Inoue (1989)]. However, the count rates for these observations were extremely low, 
and it is unclear whether or not the putative variability detections described by Ebisawa, 
Mitsuda, k. Inoue (1989) were below the effective Poisson noise limit. More recently, 
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Cui ct al. (1997) has found evidence (based upon RXTE observations) of a 3 — 12 Hz 
feature in the soft state of Cygnus X-l. We note, however, that this feature is very broad 
and is difficult to associate with a discrete feature distinct from the observed broad-band 
variability. 

Prior to the recent RXTE observations, a fairly clean division of QPO frequencies 
occurs if one ignores the 0.08 Hz QPO attributed to LMC X-l. Low frequency QPO 
(< 1 Hz) are seen in low/hard states, while high frequency QPO (~ 6 Hz) are seen in 
very high/soft states. Again, two questions arise. First, are these QPO associated with 
modes distinct from the broad-band variability? Second, if the QPO are distinct modes, 
are they associated with an accretion disc or with another component (such as a corona) 
in the BHC system? In the next subsection, we discuss a set of observed features in two 
BHC that may indeed be attributable to an accretion disc. 

3.2. QPO in 'Microquasars' 

Recently, two rather unusual and dramatic X-ray transient galactic BHC have been 
observed: GRS 1915+105 [Mirabel & Rodriguez (1994), Morgan, Remillard, & Greiner 
(1997), Chen, Swank, & Taam (1997), Taam, Chen, & Swank (1997)] and GRO J1655- 
40 [Hjcllming & Rupen (1995), Remillard ct al. (1997)]. These sources were unusual in 
that they showed powerful, highly relativistic radio jets [Mirabel & Rodriguez (1994), 
Hjellming & Rupen (1995)]; hence, they have been dubbed 'microquasars'. Furthermore, 
both have shown dramatic X-ray variability. 

As observed by RXTE, GRS 1915+105 has shown X-ray count rates up to 10 5 cts/sec, 
which, given a distance of ^ 12 kpc [Mirabel & Rodriguez (1994)], indicates a peak 
luminosity of over 10 40 erg/s [Morgan, Remillard, & Greiner (1997), Chen, Taam, & 
Swank (1997)]. As discussed in Morgan, Remillard, & Greiner (1997), GRS 1915+105 has 
also shown intense variability patterns, with the luminosity changing by factors of several 
in only a few seconds. Amidst this spectral behavior, a host of different QPO features 
have been observed, ranging in frequencies from ~ 0.1 — 10 Hz [Morgan, Remillard, & 
Greiner (1987), Chen, Swank, & Taam (1997), Taam, Chen, & Swank (1997)]. Many of 
these QPOs have multiple harmonics, and both correlations and anti-correlations with 
source luminosity have been observed. 

However, among the various QPO features seen, a high frequency QPO at 67 Hz 
has stood out because it apparently does not appreciably vary in frequency [Morgan, 
Remillard, & Greiner (1997)]. During the first epoch that this QPO feature was observed, 
it was relatively narrow (Q ~ 20), weak (rms variability ~ 0.3 — 1.6%), and varied in 
frequency by < 3%, despite factors of ~ 2 variations in source luminosity. When viewed 
in restricted bandpasses the rms variability was seen to increase with energy, with a 
maximum rms variability ~ 6% in the highest energy bandpass (~ 10 — 20 keV). During 
subsequent observations similar features were observed, always at 67+2 Hz [even on 
occasions when no soft, excess flux above a power law tail was observed in the spectrum; 
Morgan 1997, Private Communication]. 

A high frequency QPO feature has also been observed in GRO J1655-40, with frequency 
300 Hz [Remillard, et al. (1997)]. However, this QPO was only detected once, and then 
only by summing several, individual observations. It was therefore difficult to place limits 
on the stability of this feature's frequency. Like the 67 Hz feature in GRS 1915+105, 
this feature was weak (rms variability ~ 0.8%). Furthermore, it was detectable only 
in the "hardest" spectra (~ 10-20 keV). The 300 Hz feature in GRO J 1655-40 was 
somewhat broader than the 67 Hz feature in GRS 1915+105. It is also notable in that 
a very accurate mass determination, 7 ± 0.2 M Q , has been made for the compact object 
in GRO J1655-40 [Orosz & Bailyn (1997)]. 
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As we will discuss in a number of properties displayed by these features, especially 
the 67 Hz feature in GRS 1915+105, are what one expects for stable, global oscillations 
of accretion discs. We will concentrate on theories of the oscillation of the innermost disc 
regions as potential explanations of these observations. 



3.3. QPO in Active Galactic Nuclei 

As discussed in §|], the characteristic time scales are one hundred thousand to one hundred 
million times longer in Active Galactic Nuclei (AGN) as compared to galactic BHC. Thus, 
if a 0.1 s QPO period is typical of a galactic black hole, then we might expect a 3 hour to 
year-long period to be typical for AGN. Such long time scales are very difficult for current 
missions to detect, as the shorter period is on the order of the total duration of a typical 
observation, and the longer period requires many dedicated observations. However, one 
AGN, IRAS18325-5926, has shown strong evidence for a 5.5 x 10 4 s periodicity that 
was observed for 9 cycles Iwasawa, et al. (1998) |. Over these few cycles, there was no 
evidence for any strong change in phase, and the amplitude of the modulations was on 
the order of 10%. This is a somewhat stronger amplitude than observed for the 67 Hz 
and and 300 Hz features described above, and therefore are difficult to explain with the 
models described in §[s|. Further (longer) observations are required to determine to what 
extent this AGN QPO is or is not similar to the QPO seen in the microquasars. 



4. Characteristic Accretion Disc Frequencies and Time Scales 

The fact that the 67 Hz QPO feature seen in GRS 1915+105 does not appear to vary in 
frequency as a function of luminosity suggests that this feature is tied to a fundamental 
frequency or time scale in the BHC system. As this frequency is apparently independent 
of accretion rate (i.e. luminosity), the relevant time scales are likely to be gravitational 
ones. The two most important parameters for the gravitational time scales of BHC are 
the black hole mass, M, and the angular momentum, J. We shall usually normalize the 
black hole mass to a solar mass. We will characterize the black hole angular momentum 
by the dimensionless parameter a = cJ/GAl 2 , where c is the speed of light and G is 
the gravitational constant. For a Schwarzchild black hole, a = 0, while for a Kerr hole 
a < 1. Below we discuss the frequencies most relevant to orbits near the equatorial plane 
of a (possibly spinning) black hole. These frequencies are those most relevant to a thin 
accretion disc-j} 

In what follows, we shall set c = G = 1. Furthermore, physical length scales will be 
normalized to GM/c 2 (w 15 km for M = 10 Af Q ), and frequencies will be normalized to 
c 3 /GM (w 3.2 x 10 3 Hz for M — 10 Mq). In general, characteristic gravitational length 
scales increase linearly with M, whereas characteristic time scales decrease linearly with 
M. Therefore, whereas 1 ms may be a relevant time scale for a galactic black hole, days 
or months might be the relevant time scale for a massive black hole at the center of an 
AGN. There will be four main frequencies of concern to us, each a function of radius r: 
the Keplerian (i.e. orbital) frequency, the radial epicyclic frequency, the vertical epicyclic 
frequency, and the Lense-Thirring precession frequency. We describe each of these in turn 
below. 

The Keplerian frequency is the frequency with which a free particle azimuthally orbits 

f We define the equatorial plane to be the plane perpendicular to the black hole angular 
momentum axis that passes through the center of the black hole. We shall consider a cylindrical 
coordinate system (r, (j>, z) set up on this plane, with z = being in the plane. 
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the black hole. This frequency, to a viewer observing at infinity, is given by: 

( 



n = 



r 3/2 



(4.2) 



For a 10 M Q Schwarzchild black hole, this is approximately 220 Hz at r = 6. 

The radial epicyclic frequency is the frequency at which a free particle oscillates about 
its original circular orbit if it is given a radial perturbation. In classical mechanics, the 
square of the radial epicyclic frequency is given by k 2 = 4t l 2 + rdVL 2 /dr, and is exactl y 



equal to ft 2 for Keplerian orbits in an r _1 potential^] [cf. Binney fe Trcmaine (1987) |. 
In General Relativity, there is a minimum radius for which an orbit is stable to radial 
perturbations. The innermost marginally stable orbit is at r = 6 for a = 0, and moves 
inward as a increases. The epicyclic frequency is given by 

3a^ 



K 2 



9 2 



1 



G 



8a 

^3/2 



(4.3) 



and is zero at the marginally stable orbit. For a Schwarzchild black hole, k reaches a 
maximum at r — 8 and is approximately 71 Hz for M = 10 M Q . Note that for a fixed 
M, the maximum k for a ~ 1 is approximately equal to f2 at the marginally stable orbit 
for a ~ 0. This coincidence of time scales has led to alternative suggestions that the 300 
Hz feature in GRO J1655-40 can be attributed to either the Keplerian frequency at the 
marginally stable orbit or to the maximum epicyclic frequency, depending upon whether 
one adopts a ~ or o ~ 1, respectively, for this source. 

Just as the radial epicyclic frequency differs from the Keplerian orbital frequency in 
General Relativity, so does the vertical epicyclic frequency for a spinning black hole 
[ Kato (1990] , |Kato (1993)1 , |Perez, et al. (1997| . The vertical epicyclic frequency, 
is the frequency at which a free particle oscillates about its original circular orbit if it 
is given a vertical perturbation. It is equal to the Keplerian orbital frequency for both 
Newtonian gravity and Schwarzchild black holes. In general, for non-zero black hole 
angular momentum it is given by: 



n\ = 



9 2 



1 - 



4a 



3a^ 



(4.4) 



As we will discuss in § 5.2.3 , this frequency is releva nt in determining the frequ encies 
of glob al "corrugati on modes" in accretion discs [cf. Kato (1990) , Kato (1993) , Perez 
(1993), gpser (1996J . 



The final gravit ational frequency that we need to consider is the Lense-Thirring pre- 
cession frequency [ Lense fe Thirring (1918) ]. If an azimuthally orbiting ring of matter is 
tilted out of the equatorial plane, it will begin to precess due to frame dragging effects 
if the black hole has non-zero angular momentum. This precession frequency, ^lt, is 
given by the expression: 

Q LT = ^ • (4.5) 

Note that for a ~ 1 it is on the order of fl, for radii near the marginally stable orbit. 
Furthermore, Q,lt has a strong radial dependence (cf. §5.2.6). 

In Figure || we plot these various frequencies for three angular momentum parameters, 
a = 0, a = 0.5, a = 0.998. In the inner regions of accretion discs around black holes, 
these frequencies typically range from ~ 100— 1000 Hz, for M = 10 M©, with the highest 
frequencies being achieved for a ~ 1. 

f This is just the statement that a radial perturbation sends a circular orbit into an elliptical 
one, with an identical orbital period, for Newtonian free-particle orbits in an r _1 potential. 
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Figure 3. Characteristic frequencies in an accretion disc around a black hole, normalized to 
c 3 /GM « 3.2 x 10 3 Hz for M = 10 M©. The angular momentum parameter, a = cJ/GM 2 , 
where J is the angular momentum of the hole. Solid line is the Keplerian frequency, Q, dotted 
line is the radial epicyclic frequency, k, short dashed line is the vertical epicyclic frequency, Q± , 
and the long dashed line is the Lense-Thirring precession frequency, Qlt- 



Note that, at least for the Schwarzchild case, the Keplerian and epicyclic frequen- 
cies can be reasonably approximated in Newtonian gravity via the use of a "pseudo- 
Newtonian" potential. Two such potentials that have been used in the literature are 
given by: 



NW 



1 - I 

r 



12 



pw = 



2 

1 

r 



(4.6) 



[cf. |Nowak fc Wagoner (1991~)| , pazczynski fc Witta (1980)| . Both of these potentials 
have k = at r = 6, and the <&nw potential has been used for Newtonian calculations of 



modes in relativistic discs a bout Schwarzchild black holes [cf. Nowak fc Wagoner (1991 



Nowak fc Wagoner (1992)| , [Nowak fc Wagoner (1993)1 , |Nowak, et al. (1997) 



All of the above frequencies have strong radial dependences, so it is unlikely that 
any one of them can lead to a narrow feature at a discrete frequency without some 
mechanism that couples a range of radii[|. Such coupling is most likely achieved via 
hydrodynamic effects, i.e. pressure, viscosity, etc. To this end, the most relevant hydro- 
dynamical frequencies are given by the frequency of sound (for a given wavelength), oj cs , 
and the vertical Brunt- Vaisala, (i.e. buoyancy) frequency, N z . In terms of the pressure, 
P, the density, p, and the radial and vertical perturbation wavelengths, A r and X z , these 
frequencies are given by: 



2, 2 



(2nf ( A; 2 + K 2 ) 



Ki = 



oz 



(7,P)->^ 
oz 



dP 



(4.7) 



(4.8) 



For a thin disc, the sound speed is of order hil, where h is the disc half-thickness [cf. 
Bhakura fc Sunyaev (1973) | . For radial wavelengths of 0(r), the frequency of sound 
waves, ~ (h/\ r )Q, is typically much less than the Keplerian frequency. The buoyancy 
frequency, as it deals with gradients on scales of O(h), can be comparable to the ver- 
tical epicyclic frequency [cf. [Kato (1993) | . (Due to symmetry, however, the buoyancy 



f Unless, of course, there is a "special" radius, such as the marginally stable orbit radius, 
which picks out a discrete frequency. However, it then may become very difficult to achieve a 
large luminosity modulation. 
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frequency vanishes at the disc mid-plane.) Both of these frequencies play a part in the 
theories discussed in §5.2. 



5. Theoretical Models of Disc Oscillations 

5.1. Historical Review 

The study of oscillations relevant to accretion discs around black holes has a history 
that dates back more than 30 years. It was quickly realized that the a— disc model of 
jhakura fc Sunyaev (1973)| shows both viscous and thermal instabilities for high lumi- 
nosities above 0(10% LEdd) jhakura fe Sunyaev (1976) |. These instabilities have often 
been invoked to explain the intense, broad band variability seen in BHC low/hard states^. 

As discussed in § |3.1[ until recently there has been little observational evidence for 
discrete (possibly stable) periodic modes in accretion discs about black holes. However, as 
was first realized by Kato fc Fukuc (1980)| , not only can accretion discs support discrete 
oscillatory modes, but also the effects of General Relativity modify the mode spectrum 
and determine the regions in which modes can be trapped. By considering adiabatic 
perturbations of an isothermal disc, Kato & Fukue (1980) showed that the roll-over of 



k to zero at the marginally stable orbit leads to a cavity that effectively traps acoustic 



.i 



modes (i.e. p— modes) at the disc inner edge (cf. 

Later, Okazaki, et al. (1987)| considered isothermal perturbations of isothermal discs 



(i.e. 



and showed that it was also possible to trap what are essentially internal gravity modes 
g— modes) near the epicyclic frequency maximum (cf. § |5.2.2 ). Vishniac & Diamond 
i considered travelling wave versions of g— modes that had an azimuthal dependence 
oc exp[?m0], with m = 1. These modes were invoked as a possible mechanism for 
angular momentum transport in an accretion disc. [Note, however, that the strictly 
Newtonian modes considered by Vishniac & Diamond (1989) were not trapped modes, 
and furthermore they had very small vertical wavelengths, i.e. X z <C h.] 

In a series of papers Nowak & Wagoner (1991,1992,1993) adopted a Lagrangian pertur- 
bation approach [cf. Lynden-Bell & Ostriker (1967), Freidman & S chutz (1978a, b)] to the 



stud y of disc oscillations. T hey showed that the acoustic modes of Kato fc Fukue (1980) 
and Okazaki, et al. (1987) are high and low frequency limits, respectively, of the same 



dispersion relationship, which itself is the strong-rotation limit of the simplest form of 
the "helioseismology" dispersion relationship [cf. Hincs (1960) [. Hence, they gave these 
mod es t he name "diskoseismology" . These works used a pseudo-Newtonian potential 
(eq. L6) to mimic the effects of General Relativity; however, they also generalized the 
calculations to adiabatic perturbations of non-isothermal discs. As discussed by Frei- 
dman & Schutz (1978a, b), by adopting a Lagrangian perturbation approach one can 
calculate a conserved (in the absence of dissipative forces) "canonical energy" for the 
modes||. Accurately determining the canonical energy allows one to determine the effects 
of various dissipative effects, such as viscosity [ |Nowak fc Wagoner (1992) |, and allows 



one to consider various excitation mechanisms for the mode ]Nowak fc Wagoner (1993) 
Nowak, et al. (1997)|. Nowak & Wagoner (1993,1997) also considered ways in which the 



modes might lead to actual luminosity modulations (cf. § |5.2.5] ). 

Ipser & Lindblom (1992) developed a scalar potential formalism for calculating modes 

f It should be noted, however, that as shown in Fig. |l| and discussed in 
Nowak f!995h low/hard state luminosities are nominally below the instability limits of 
jhakura fc Sunyaev (1976) . 



X As Freidman & Schutz (1978a,b) point out, assigning an energy to a mode in a rotating 
medium is an extremely subtle issue, and many prior works were incorrect in determining this 
quantity. 
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of rotating systems in full General Relativity. Gradients of the scalar potential are 
related to the Lagrangian perturbation vector in both the pseudo-Newtonian formalism 
[Nowak & Wagoner (1992)] and the fully General Relativistic formalism [Perez, et al. 
(1997)]. Perez, et al. (1997) used this formalism to consider the fully relativistic version 
of g— modes. In addition, a third class of modes, the c— mode s (cf. jj|5.2.3 ) recently have 



been identified and described with fully relativistic formalisms [Kato (1990), Kato (1993) 
Perez (1993), Ipser (1996).] 



5.2. 'Diskoseismology' 

The main distinguishing feature of most of the theoretical modes described in the works 
cited above is that the mode frequencies predominantly depend upon fundamental grav- 
itational frequencies and are not strongly effected by hydrodynamic processes if h <C r. 
That is, for a thin disc, changes in density and/or pressure (which we assume to be 
correlated with luminosity) may effect the physical extent of a mode or its amplitude, 
but they do not greatly effect the mode's frequency. Furthermore, the expected mode 
frequencies, ^ 100 Hz, are comparable to the gravitational frequencies in the innermost 
regions of the accretion disc. 

From these standpoints, the QPO features discussed in §[0] are not good candidates 
for 'diskoseismic' modes. They are of very low frequency, and furthermore several, such 
as the 3 — 8 Hz QPO in Nova Muscae [Miyamoto, et al. (1994)], have been seen to 
vary in frequency by fractionally large amounts. Likewise, the low frequency (~ 0.1 — 10 
Hz) QPO seen in both GRS 1915+105 and GRO J1655-40 have been observed to be 
very highly variable in frequency [Morgan, Remillard, & Greiner (1997), Remillard, et 
al. (1997)]. 

However, the 67 Hz feature observed in GRS 1915+105 has been consistently observed 
at 67+2 Hz, despite the fact that this source's luminosity has varied by factors of two or 
more over the epochs during which this feature has been observed. Although there has 
been only one detection of the 300 Hz feature in GRO J1655-40, it was consistent with 
being at a steady frequency. The 300 Hz, if attributable to a disc time scale, indicates a 
phenomenon occurring very close to the disc inner edge. For these reasons, we identify 
both of these features as candidate 'diskoseismic' modes. 

For pseudo-Newtonian potential calculations of diskoseismic modes, one defines a La- 
grangian perturbation vector that describes the displacement of a fluid vector relative 
to its unperturbed path [cf. Freidman & Schutz (1978a)]. All perturbation quantities 
compare the displaced, perturbed fluid element to the (moving) fluid element in the 
unperturbed flow. If we take £ to be the displacement vector, then the Lagrangian varia- 
tion of the velocity, Av, is given by d^/dt [Freidman & Schutz (1978a)]. [As in Nowak & 
Wagoner (1991), we use A to denote a Lagrangian perturbation.] The Lagrangian per- 
turbation of the density is derived from mass conservation, and is given by Ap = — pV 
For adiabatic oscillations, Lagrangian perturbations of the pressure are related to La- 
grangian variations of the density by AP/ P = ^Ap/ p, where 7 is the adiabatic index 
[cf. Freidman & Schutz (1978a), Nowak & Wagoner (1991)]. 

Calculationally, it is somewhat easier to determine the mode structure by defining a 
scalar potential [Ipser & Lindblom (1992), Nowak & Wagoner (1992)], 6V(r,z) = SP/p, 
where here 5 refers to an Eulerian perturbation||. The potential SV is seen to be the 
Eulerian perturbation of the enthalpy, and furthermore the Lagrangian perturbation is 
related to gradients of this potential [Nowak & Wagoner (1992)]. If one looks for modes 

f An Eulerian perturbation is one that compares the perturbed fluid to the unperturbed fluid 
at a, fixed coordinate [Freidman & Schutz (1978a)]. 
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oc exp[— i(m4> + at)], where t is time and a is the inertial mode frequency, then in a WKB 
approximation the modes satisfy a dispersion relationship of the form 

[u 2 - T(r)n 2 ] [u 2 - n 2 ) = u; 2 c 2 s k 2 , (5.9) 

where lu = a + mil is the corotating frequency, T(r) is a slowly varying separation func- 
tion (which comes from separating the basic fluid perturbation equation into radial and 
vertical components), £1 and n are the Keplerian rotation frequency and epicyclic fre- 
quency, respectively, c s is the speed of sound, and k r = 2ir /X r is the radial wavenumbcr. 

This dispersion relation is essentially the same as the simplest form of the helioseis- 
mology dispersion relation, except that k takes the role of the buoyancy frequency, and 
^/T(r)0 takes the role of the 'acoustic cutoff' frequency [cf. Hines (I960)]. Nowak & 
Wagoner (1991, 1992) showed that there are two general classes of solutions to this dis- 
persion relation: high and low frequency. We identify the high frequency (lu 2 ^ k 2 ) 
modes with acoustic p— modes [cf. Kato (1980), Nowak (1991), Perez (1993)], and the 
low frequency (lu 2 ^ k 2 ) modes with internal gravity g— modes [cf. Okazaki et al. (1987), 
Vishniac & Diamond (1989), Nowak & Wagoner (1992), Perez, et al. (1997)]. Very qual- 
itatively, gravitational effects set certain fundamental oscillation frequencies, which are 
then modified by pressure effects. If the pressure forces 'assist' the gravitational forces, 
then a high frequency is achieved. If the pressure forces 'retard' the gravitational forces, 
then a low frequency is achieved. (We will return to this qualitative notion below when 



we consider mode excitation, cf. §5.2.4 



The fully General Relativistic version of the diskoseismology equations leads to a WKB 
dispersion relationship of a very similar form. As was shown in Perez, et al. (1997), for 
the relativistic equations a WKB analysis still allows approximate separation of the 
governing equation into radial and vertical dependences. The radial component of the 
fluid perturbations satisfies the WKB relation 



d 2 W 
dr 2 



M'l — i - I 



2 )W = 0, (5.10) 



where a{r) is inversely proportional to the speed of sound at the mid-plane and ^(r) is 
a slowly-varying separation function. The eigenfunction W(r) is proportional to a radial 
derivative of the potential SV (again defined as SP/p) and to the radial component of 
the Lagrangian fluid displacement [Perez, et al. (1997)]. 

From the radial WKB equation ( 5.10| ) one can identify three classes of modes that 



are trapped. The p— modes, defined by ty(fl±/uj) < 1, are trapped where or > n 
g— modes, defined by ^(il±/ui) 2 > 1, are trapped where lu 2 < n 2 . The third class 
of modes (unique to the fully relativistic treatment), the c— modes, are defined by 

(f2 j_ / 'lu) 2 = 1 and may in principle be trapped in either region. [The separated equa- 
tion governing the vertical component involves a more complicated second-order linear 
operator which contains the vertical buoyancy (Brunt- Valsala) frequency, N z .] 

The frequencies of all classes of modes are proportional to 1/M, but their dependences 
on the angular momentum of the black hole are quite different. In principle this would 
allow one to measure the angular momentum of the black hole if more than one type of 
mode were to be detected in the same source. Alternatively, if one could infer the mass 
of the black hole through the motion of a companion which feeds the disc, such as has 
been done for GRO J1655-40 [Orosz & Bailyn (1997)], one could determine the angular 
momentum of the black hole by a single mode observation (given a correct identification 
for the class of mode observed). Below we briefly describe each class of mode, and we 
present representative examples in Figures ^|-^[ 
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Figure 4. A representative m — p— mode. Here we show the Eulerian potential, 5V, as well 
as the radial component of the Lagrangian perturbation vector as a function of radius r, and 
normalized vertical coordinate, x. [Taken from Perez (1993).] 

5.2.1. P-modes 

As discussed by Kato (1980) and Nowak & Wagoner (1991), p— modes are trapped in 
regions where oj 2 > k 2 . For m = 0, this leads to a narrow trapping region (Ar mo de "C 1) 
at the disc inner edge where k 2 goes to zcrofj]. These modes have frequencies a factor 
of a few less than the radial epicyclic maximum, and tend to have £ r <; £ z (i.e. ra- 
dial perturbations greater than the vertical perturbations) [Nowak & Wagoner (1991)]. 
Due to their narrow confinement very near the disc inner edge [where the 'no-torque' 
disc boundary condition leads to very little luminosity modulation; Shakura & Sunyaev 
(1973)], we do not expect these modes to have significant luminosity modulation. We 
present an example of a p— mode in Fig. ^. 

5.2.2. G-modes 

Internal (gravity) modes are trapped where uj 2 < k 2 , in the region where k achieves 
its maximum value and where the disc is hottest. The lowest modes can have significant 
vertical displacements (£ z <; £ r ) and relatively large radial extents Ar « GM/c 2 . We 
believe that these modes will produce the greatest luminosity modulations in the disc; 

f Trapped p— modes might also exist in the large, outer region of the disc; Silbergleit 1997, 
Private Communication. 
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thus, they may be the most observable class of modes [Nowak & Wagoner (1992), Perez, 
et al. (1997)]. 

To analytically approximate the eigenfunctions and eigenfrequencies of the lowest 
g— modes, WKB solutions to the separated radial and vertical equations of fluid per- 
turbations can be obtained. From the symmetry of the governing equations, it is suffi- 
cient to consider eigenfrequencies a < and axial mode integers m > 0. The resulting 
frequencies, / = —a/2ir 1 of the lowest radial (m = 0) g— modes are given by 



/ = 714 (1 - e nj )(M Q /M) F(a) Hz, 

n + i\ h 



j + 5 J r 



(5.11) 



Here F(a) is a known, monotonically increasing function of the black hole angular mo- 
mentum parameter a, with F(a = 0) = 1 and F(a = 0.998) = 3.44. The properties 
of the disc enter only through the small correction term e n j, which involves the disc 
thickness 2h(r) and the radial (n) and vertical (j) mode numbers, with S ~ 1. Typically 
h/r ~ 0.1L/Z,Edd for a radiation-pressure dominated optically thick disc region, where 
L/LEdd is the ratio of the luminosity to the Eddington (limiting) luminosity^. [Higher 
axial g— modes with m > have a somewhat different dependence on a than the radial 
modes; Perez, et al. (1997)]. 

The 67 Hz feature observed in GRS 1915+105 has been associated with a diskoseismic 



g— mode [Nowak, et al. (1997)]. If this association is correct, then equation (5.11) 
predicts a black hole mass of 10.6 M Q , if the hole is non-rotating, up to 36.3 M Q , if 
t he ho le is m aximally rotating. (Further aspects of this identification are discussed in 



§ |5.2.4| , § p.2.5| below.) If we identify the 300-Hz feature observed in GRO J1655-40 with 
the fundamental g— mode oscillating in an accretion disc surrounding a 7.0 Mq black 
hole [as deter mined from spectra of the companion star; Orosz & Bailyn (1997)], then 



equation 5.11 implies that its angular momentum is 93% of maximum. 



5.2.3. C-modes 

To analytically obtain the approximate eigenfrequencies for the c— modes, one can solve 
the separated WKB perturbation equations in the regime 1 &(Q±/uj) 2 = 1 [Lehr, Wagoner, 
& Silbergleit (1997), Private Communication]. Here we consider c— modes that are non- 
radial (m > 1), nearly incompressible oscillations trapped in the very innermost region 
of the disc with eigenfrequencies |er| ~ mQ,(r c ) — Q±(r c ), where r c is the upper radial 
bound of the mode. For m = 1 and [a/r'l <C 1 the fundamental c— mode eigenfrequency 
is approximately the Lense-Thirring frequency evaluated at the radius r c : 

H^- (5-12) 

' c 

Note that for the m = 1 mode, r c is typically no greater than 10% larger than the 
marginally stable orbit radius. The physical structure of this mode resembles a tilted 
inner disc which slowly precesses about the black- hole spin axis (cf . Fig ||) . Since the 
mode is nearly incompressible, there is little temperature or pressure fluctuation in the 
disc, and hence a c— mode may produce very little intrinsic luminosity modulation (cf. 



5.2.5). The mode may be observable, however, since the projected area of the inner 



f Although the frequencies of the modes do not have a strong dependence upon L/L^dd, 
the modes are no longer effectively trapped at large radii if L/L^dd 0.3. For such high 
luminosities, the outer evanescent region for the g— modes becomes narrow in radius, thereby 
allowing g— modes to effectively couple to travelling waves in the outer regions of the disc. The 
modes thus "leak" energy to larger radius. 
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Figure 5. A cross section of a disc with an m = g— mode. The disc thickness, as a function 
of radius, has been normalized to 1. Shaded contours correspond to the Eulerian pressure 
perturbations normalized to the local unperturbed pressure values. (The apparent asymmetry 
about the disc mid-plane is a numerical artifact, due to inadequate resolution of the contouring 
routine.) 

disc changes with time. This allows the possibility of coronal emission being Compton 
reflected and modulated in the disc inner regions. Calculations to determine the extent 
of the resulting modulation are currently being undertaken [Nowak & Reynolds (1997), 
Private Communication] . 

Note that GRO J1655-40 does have a high inclination, ~ 70°, to our line of sight 
[Orosz & Bailyn (1997)] making this a potentially promising mechanism for producing 
the observed 300 Hz feature. (Also, Compton reflection features would peak near ~ 30 
keV, which is also consistent with the observations.) If we associate the 300 Hz feature 
with a c— mode in a disc about a 7 M black hole, then a « 0.8 for this source. Note 
that this is different from both the prediction made from assuming a g— mode in this 
source, as well as the prediction of Cui, et al. (1998) who associated this feature with 
Lense-Thirring precession at a single radius (cf. |5.2.6| ). In Fig. ^ we show the c— mode 
frequency as a function of the angular momentum parameter a. 

5.2.4. Excitation and Damping Mechanisms 

It is possible to use a parameterized stress tensor to estimate the effects of turbulent 
viscosity on the g— modes [Nowak & Wagoner (1992,1993)]. The canonical energy of 
a radial mode is E c ~ o- 2 p(£ 2 + dV, where here dV is the volume occupied by 
the mode. Isotropic turbulence produces a rate of change dE c /dt = —E c /t, with r ~ 
| aa [h 2 /\ 2 + h 2 /Xl] Here A r and A z , respectively, are the radial and vertical mode 
wavelengths, and a is the standard Shakura-Sunyaev a— parameter [Shaukura & Sunyaev 
(1973)]. The corresponding quality factor is given by 

Qjn- 1 = (\o-\t)- 1 ~ [j 2 + (h/r)n 2 ] a, (5.13) 
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Figure 6. Left: Schematic drawing of a c— mode in an accretion disc (arbitrary amplitude). 
Right: Frequencies of the fundamental m = 1 c— mode frequency as a function of the black hole 
angular momentum parameter a = cJ/GM 2 . Frequencies are normalized to c 3 /GM « 3.2 xlO 3 
Hz for M = 10 M Q . 



as X z ~ h/j and X r ~ vhr/n, where j and n are of order of the number of vertical and 
radial nodes in any particular eigenfunction. Thus, for a <C 1, we can have high mode 
Q [cf. Nowak & Wagoner (1993); Nowak, et al. (1997)]. 

The above estimates are for isotropic viscosity. If the turbulence does not efficiently 
couple to the vertical gradients of the modes, then the mode Q value is increased by a 
factor^ (j\ r /h) 2 [cf. Nowak & Wagoner (1993)]. Aside from damping modes, turbulence 
can also potentially excite modes. Velocity perturbations in the disc, 8v, are made up of 
a mode component, Svm, and a turbulent component, Svt- Viscous damping arises from 
terms of the form SvMifivTj , while mode excitation arises from terms of the form SvTiSvxj 
[Nowak & Wagoner (1993)]. It is possible to make simple estimates of the magnitude 
of the turbulent excitation, and balance this against the turbulent damping [Nowak & 
Wagoner (1993)]. The modes are excited to an amplitude of |£ 2 | ~ a(h/X r ) 3 ^ 2 h, and 
~ ay^Xr/h h, for isotropic and anisotropic viscosity, respectively [Nowak & Wagoner 
(1993)]. If turbulence is playing the dominant role in damping and exciting the modes, 
then we have the following constraints. For isotropic turbulence, the Q value is large 
only for a«l; however, this implies a correspondingly small amplitude. For anisotropic 
viscosity, a larger mode amplitude and a higher Q is achieved for a given a. However, 
in either case it is likely very difficult to achieve amplitudes as large as required to agree 
with the observations of GRS 1915+105 and GRO J1655-40 (cf. § fjl| ). 

Another potential excitation mechanism for g— modes is 'negative radiation' damping 
[Nowak, et al. (1997)]. As a first approximation, the g— modes are taken to be adiabatic. 
In reality, one expects there to be small entropy changes due to various effects, the most 
notable one being radiative losses. If one has a radiation pressure dominated atmosphere, 
as is likely in high-luminosity discs, one properly should use 

P P 

AP = 7— Ap + 7— As , (5.14) 
P s 

where s is the specific entropy. The effects of this term for a radiation pressure dominated 
atmosphere can be approximated by writing an 'effective' adiabatic index, 7', to be 
substituted into the relationship between the Lagrangian variations of the pressure and 
density. Nowak, et al. (1997) showed this term to be 

Y « 7(1 + i^rj^ = 7(1 + io!) , (5.15) 
where r es is the scattering depth of the disc and h is the disc half thickness. In the 
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above, our ignorance of the disc's true vertical structure is subsumed into the parameter 
a' , which one expects to be of 0(a) for a radiation pressure-dominated disc. 

It is the imaginary component of the effective adiabatic index that leads to radiative 
damping for p— modes but conversely leads to excitation of the g— modes. One can see 
this by perturbing the approximate dispersion relationship lo 2 ss k 2 ± c 2 k 2 , where the 
+ is for p— modes and the — is for g— modes. Perturbing the frequency and the sound 
speed together, we have: 



where Sc 2 ~ ia'c 2 . Thus, the imaginary component of 7' leads to damping of p— modes 
and excitation of g— modes. Qualitatively, as the p— mode compresses gas, radiation leaks 
out, thereby providing less restoring force from the compressed region, and therefore 
leads to damping. As discussed above, the g— mode is in some sense opposite as the 
pressure works against the gravity. Hence, radiative leaking leads to less effective pressure 
restoring forces and thereby to g— mode growth. 

5.2.5. Luminosity Modulation 

For most simple a-models, the energy generation rate per unit volume in the accre- 
tion disc is approximated by aP(r, z)£l(r). The local pressure is modulated by p— and 
g— modes (and to a much lesser extent by c— modes); thus, there is the possibility that 
these modes can lead to an observable luminosity modulation. Furthermore, the modes 
not only perturb the pressure, but they also perturb the locations of the disc boundaries. 
As discussed in Nowak, et al. (1997), the integrated variation of the energy generation 
rate is therefore 



where P'(r',z') = P(r 7 z) + AP(r, z), r' = r + £ r (r), z' = z + £ z (z), and 77, ro, z n are 
the disc boundaries (ro can be taken to go to 00 without loss of generality). 

The <?-modes have the most promising combination of pressure modulation and physical 
extent, which makes them the most likely candidates to be observed by direct luminosity 
modulation of this sort. However, for standard disc models, only a few percent of the 
total luminosity is generated in the region of the disc where the g— modes can exist. Thus, 
we expect only of 0(1%) rms modulation due to the g— modes. This is consistent with 
the observations of GRS 1915+105 and GRO J1655-40; however, as shown in Nowak, 
et al. (1997), even this modest luminosity modulation requires large mode amplitudes. 
Specifically, £ z ~ h is required^. 

Qualitatively, as the g— modes exist in the inner, hotter disc regions, one expects the 
modulated luminosity to be harder than the average disc luminosity. However, both the 
67 Hz and 300 Hz QPO seem to be strongest in the power law tail, which is problematic 
for the g— mode interpretation of these modes if one associates the power law tail with a 
corona distinct from the accretion disc. 

On the other hand, although one expects the c— modes to have relatively little luminos- 
ity modulation due to pressure fluctuations, these modes may be effective at modulating 

f Note, however, that Nowak, et al. (1997) only considered Schwarzchild black holes. As a 
disc about a Kerr black hole releases a fractionally greater amount of its energy in the g— mode 
region, one may only require a mode amplitude a factor of several smaller. 




(5.16) 





(5.17) 
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reflected emission. If the power law tail is due to a corona above a cold accretion disc, 
then part of this power law will be Compton reflected [Magdiarz & Zdziarski (1995)], 
and hence may be modulated by the warping associated with the c— modes. Further- 
more, one expects the reflected emission to peak in the 20 — 30 keV range [Magdiarz & 
Zdziarski (1995)], which is consistent with the observations of both GRS 1915+105 and 
GRO J1655-40. 

5.2.6. Alternative Models 

The diskoseismology g— mode interpretation of the 67 Hz and 300 Hz has several points 
in its favor. First, it reproduces the correct time scales. Second, the g— mode frequencies, 
like the observed 67 Hz frequency, are expected to be relatively insensitive to luminosity 
fluctuations. Third, mode Q of O(20) is achievable with a reasonable a— parameter. 
Fourth, there is an identified mechanism (negative radiation damping) for exciting the 
modes. Fifth, the <?— modes are capable of modulating the luminosity of the disc, possibly 
as much as the observed ~ 1%. Furthermore, we qualitatively expect the modes to be 
harder than the mean disc spectrum. 

The are two major arguments against the g— mode interpretation. First, one requires 
a large mode amplitude (although this might be reduced slightly if one considers near 
maximal Kerr holes). Second, the observed QPO features appear strongest in the power 
law tails, which are not necessarily associated with a disc component. Aside from the 
c— modes discussed above, there are several other alternative hypotheses to the g— modes. 

One possibility is fluctuations in the location of the sonic point [Honma, Matsumoto, 
& Kato (1992)], which can lie just inside of the marginally stable orbit, where disc orbits 
are going from circular orbits to near free-fall trajectories. As discussed by Honma, 
Matsumoto, & Kato (1992), the time scales associated with this 'pulsational instability' 
of the sonic point are fairly rapid. Again, one needs to show that sufficient luminosity 
modulation is possible with such a mechanism. 

Milsom & Taam (1997) have numerically shown that for some accretion disc mod- 
els, one finds acoustic p— modes that, although not trapped, peak at the location of the 
epicyclic frequency maximum^. Furthermore, their frequency is that of the epicyclic fre- 
quency maximum. Many predictions of the p— mode hypothesis, such as the mass and 
angular momentum of the central object, are identical to those made by the g— mode 
hypothesis. As for the g— modes, it is also difficult to achieve (9(1%) luminosity modula- 
tion, and likewise it may also be difficult to produce features as hard as those observed 
in GRS 1915+105 and GRO J1655-40. 

More recently, Cui, Zhang, & Chen (1998) have proposed that the observed QPO fea- 
tures are associated with Lense-Thirring precession. Specifically, these authors choose the 
QPO frequency to be the Lense-Thirring frequency at the radius at which the accretion 
disc effective temperature is maximum, and thus derive a ~ 1 for both GRS 1915+105 
and GRO J1655-40. There are a number of objections to this picture. First, choosing the 
radius to be that of the effective temperature maximum is somewhat arbitrary, and one 
might expect that the radius at which r 2 .F 7 is maximized, where F 7 is the disc photon flux 
over the observed energy bands, is a more natural choice. Second, neither F 7 nor r 2 F-y 
are sharply peaked functions, therefore it is extremely unlikely that one can generate a 
feature with Q ~ 20, as was observed for th e 67 H z feature in GRS 1915+105. [From this 



point of view, the c— modes discussed in £5.2.3 are the natural, 'global' features to be 



t Qualitatively, one can understand this by considering the p— mode dispersion relationship, 
which is (uj 2 — k 2 ) ~ c 2 k 2 . The p— mode group velocity is seen to go to zero at uj 2 = k 2 . 
Thus p— modes with the maximum epicyclic frequency "stall" at the location of the epicyclic 
frequency maximum. 
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associated with Lense-Thirring precession.] Third, due to frame dragging effects in the 
presence of viscous forces, one expects that the inner region of the disc will be flattened 
and constrained to the equatorial plane on a time scale of 0(ls) [Bardeen & Petterson 
(1975)f§ Fourth, Cui, Zhang, & Chen (1998) do not identify any luminosity modula- 
tion mechanism. Finally, Cui, Zhang, & Chen (1998) do not demonstrate that there is 
a viable excitation mechanism for the modes [although they suggest that the radiative 
warping mechanism of Pringle (1996) may be at work]. However, as discussed above, the 
c— modes, which are qualitatively similar to the Lense-Thirring precession suggested by 
Cui, Zhang, & Chen (1998), may yet provide a viable explanation for the observed high 
frequency QPO features. 



6. Summary 

As was discussed in §5.1, the study of stable oscillations in accretion discs has a long 
and rich history. However, it was not until the advent of Ginga in the late 1980's and of 
RXTE in only the past few years that this field has also become an observational one. 

Prior to the launch of RXTE, there were only a few observations of QPO features in 
BHC. Most of these features were of low frequency ^ 10 Hz (this was likely mainly due 
to instrumental limitations), were not seen during more than one epoch, and were often 
broad (Q 10) and often variable in frequency. Few theories have been put forth that 
adequately describe these observations. 

In the past two years since the launch of RXTE, a wealth of new observational infor- 
mation has been obtained. Again, a number of low-frequency, often broad and highly 
variable features have been observed. Among the wealth of features seen with RXTE, 
two high frequency features stand out: the 67 Hz feature in GRS 1915+105 and the 300 
Hz feature in GRO J 1655-40. The former is notable for its steady frequency, whereas the 
latter is notable for its high frequency which suggests very strongly that it comes from 
very close to the probable 7 M© black hole in this system. 

The high frequency of both features and the stability of the 67 Hz feature suggests 
that these QPO might be related to stable oscillations in the inner regions of accretion 
discs. We have described a class of theories, which we refer to as 'diskoseismology', that 
might offer an explanation for these observations. The main motivations for attributing 
these features to diskoseismic modes are that these modes: 1) are related to a "natural" 
frequency in the disc (i.e. the maximum epicyclic frequency); 2) their spectra are ex- 
pected to be characteristic of the inner, hottest regions of the disc; 3) their frequencies 
are relatively insensitive to changes in luminosity; and 4) they have low rms variability. 

This latter feature, although in agreement with the observations, is the strongest con- 
straint. These modes cannot be applied to systems that show <; 10% rms variability 
over a wide range of energy bands. Furthermore, it may be difficult to explain the ob- 
served spectral hardness of the QPO. To these ends, one needs to begin to consider 
more detailed dynamical disc models that address the production of the hard radiation. 
Also, g— modes have been the major focus of recent study; however, c— modes, which are 
related to Lense-Thirring precession of warped accretion discs, may offer a better expla- 

X Cui, Zhang, & Chen (1998) mistakeniy ciaimed that this time scaie was very long; however, 
they were quoting the time scaie for the black hole's angular momentum to align with the binary 
orbital angular momentum (which is of order millions of years). The time scale for the inner 
region of the disc to flatten into the orbital plane is significantly shorter. Note that the c— modes 
should also be susceptible to damping by the 'Bardeen-Petterson' effect; however, as they are 
global modes we expect the damping time scale to be slightly longer than that suggested by 
Bardeen & Petterson (1975). 
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nation of some of the data. Finally, in all of the above work magnetic fields have been 
neglected. As shown by Balbus & Hawlcy (1992), even weak magnetic fields can play an 
important dynamical role. The incorporation of magnetic fields into the diskoseismology 
calculations is one of the next major steps that needs to addressed. 

Even if the features seen in GRS 1915+105 and GRO J1655-40 do not turn out to be 
diskoseismic modes, they point out two important lessons. First, BHC systems can pro- 
duce relatively stable, high-frequency features. Second, the Rossi X-ray Timing Explorer 
is capable of detecting and characterizing these features despite their weak variability. In 
a very real sense, despite nearly twenty years of research, the study of stable oscillations 
in black hole accretion systems has just begun. 
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